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It is shown analytically that there exists a natural basis in terms of which the nonperturbative 
time evolution of an important class of driven four-level systems in the strong-coupling regime 
decouples and essentially reduces to the corresponding time evolution in the weak-field regime, 
exhibiting simple Rabi oscillations between the different relevant quantum states. The predictions 
of the model are corroborated by an exact numerical calculation. 
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I. INTRODUCTION 

S . 

The dynamical behavior of quantum systems driven by external time-dependent fields has attracted considerable 
' interest in recent years due, in part, to the great variety of phenomena that have been theoretically predicted and 
£SJ , experimentally observed when the system conditions are conveniently chosen [1-3]. For instance, in the field of 
quantum optics the quantum interference effects induced by the coherent external fields can lead to phenomena such 
as coherent population trapping [4] (even in the nonperturbative regime [5]), electromagnetically induced transparency 
[6], or lasing without inversion [7]. In atomic systems an external laser field can induce interesting processes such as 
harmonic generation [8] and multiphoton excitation and ionization [9]. 

The theoretical treatment of a quantum system exposed to a strong time-dependent field requires specific nonper- 
turbative methods. A first comprehensive theoretical study of the effects of a strong oscillating field on a two-level 
quantum system was carried out by Autler and Townes [10], who making use of Floquet's theorem [11] derived a 
solution in terms of infinite continued fractions to investigate the effect of an rf field on the J = 2 — > 1 /-type doublet 
microwave absorption lines of molecules of gaseous OCS, obtaining good agreement with the experimental results. 
In another important paper Shirley [12] used also the Floquet's theorem to develop a general formalism for treat- 
<*S j' ing periodically driven quantum systems. Using this formalism, which replaces the solution of the time-dependent 
, 1 * . Schrodinger equation with the solution of a time-independent Schrodinger equation represented by an infinite matrix, 
he obtained closed expressions for time-average resonance transition probabilities of a strongly-driven two-level sys- 
tem. More recently, a variety of approaches have been proposed to deal analytically with strongly driven two-level 
systems [13-19]. Three- and four-level systems driven by intense laser fields has also been treated analytically [5,20]. 

In the numerical description of realistic multi-level atoms and molecules in intense laser fields the Floquet theory 
and, more recently, the R-matrix-Floquet approach [21], have also proved to be particularly useful. These formalisms 
have been used in studies of atomic spectroscopy [22], laser-assisted electron-atom scattering [23], harmonic generation 
[24] , periodically kicked Rydberg atoms [25] and multiphoton excitation and ionization of atoms and molecules [26-30] . 

In this work we are interested in an important class of driven four-level systems. Specifically, we consider a four-level 
system consisting of two doublets (see Fig. 1). This system has been previously studied in the context of coherent 
population transfer [31] and tunneling dynamics [20]. In the present work we will show that there exists a natural 
basis in terms of which the nonperturbative time evolution of the system in the strong-coupling regime decouples and 
essentially reduces to the corresponding time evolution in the weak-field regime. 

The splittings of the two lower states (|1), |2)) and the two upper states (|3), |4)) will be denoted as Aq and A ', 
respectively. These splittings are much smaller than the separation A between the doublets. Such a level configuration 
is commonly encountered in quantum double- well potentials, which in turn are relevant in the description of numerous 
processes in molecular and solid-state systems. For instance, this model can describe the tunneling dynamics of the 
inversion mode of the ammonia molecule [1,32], intermolecular proton transfer processes [33], or the effect of a driving 
laser field on the tunneling dynamics of low-lying electrons in quantum semiconductor heterostructures [34,35]. 

The external periodic field (of amplitude E and frequency oj) will induce transitions between states |1) <-> |2), 
3) <-> 4), |1) <-> |4), and |2) <-> |3), with corresponding coupling constants f2i2, ^34, ^14, and ^23, where f2y = E/i^ 
and fiij is the dipole matrix element between states |i) <-> The Hamiltonian of the system reads 
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FIG. 1. Energy diagram and coupling constants of the four-level system considered in the present work. 
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H =} EiOa - ^12 COS (uit) ((712 + 021 ) - ^34 COS (u)t) ((7 34 + <7 43 ) 



-^14 COS (uit) ((7i4 + O-41) - f2 2 3 COS (wi) (cr 2 3 + 0-32) , 



(1) 



where = is the transition operator, £7j is the energy of state |i) in the absence of the periodic force, and we 
take h = 1 throughout the paper. 

We shall assume that, as is usually the case, the dipole matrix elements Hij between states lying within a given 
doublet are much larger than the corresponding dipole matrix elements connecting states lying in different doublets. 
Under these circumstances, the states within a doublet are much more strongly coupled by the external field than those 
lying in different doublets, i.e., f2i2, ^34 S> f2i4, ^23- If the driving field is quasiresonant with the allowed transitions 
between the lower and upper doublets (u> « A) and weak enough so that any coupling constant Of,- is much smaller 
than the frequency ui (weak-coupling regime), the contribution of the far-off-resonant transitions |1) <-> |2) and 
1 3) <-» |4) turns out to be negligible. Under these conditions one can invoke the rotating wave approximation (RWA) 
and the dynamical evolution of the system becomes governed by a Hamiltonian which, in the rotating frame, takes 
the simple form H' = H[ 4 + with 



and a similar expression for H' 2 ^ replacing 1 — > 2 and 4^3. Thus, in the weak-field regime, the time evolution 
consists of usual Rabi oscillations between |1) <-» |4) and 2} <-> |3). As the strength of the external field increases, 
the contribution of the strong nonresonant-transitions |1) <-> \2) and |3) <-> |4) becomes increasingly important so 
that, eventually, the system enters an interesting regime in which the intra-doublet transitions become strong while 
the corresponding inter-doublet transitions remain weak. This is a nonperturbative strong-coupling regime where the 
RWA is no longer valid and one has to deal with two weak and two strong transitions. Under these circumstances 
all the states become coupled and the dynamical evolution becomes, in general, rather involved. As we will show, 
there exists, however, a natural basis in terms of which the time evolution of the system is essentially the same in 
both the perturbative and nonperturbative regimes. This basis thus provides a unified description of the weak- and 
strong-coupling regimes. 

Our approach is not directly based on Floquet theory, rather it relies on a suitable time-dependent unitary transfor- 
mation which allows the (intra-doublet) strong contributions to be conveniently absorbed into renormalized physical 
parameters. However, a connection can be established between these two approaches. Floquet states and quasiener- 
gies, which are, respectively, the eigenstates and eigenvalues of the hermitian operator (H — id/dt), become in the 
time-independent case indistinguishable from the usual stationary states and energies. Thus, if one performs, as we 
shall do, a unitary transformation U(t) to a rotating frame in which the transformed Hamiltonian becomes, to a good 
approximation, time- independent, then finding the Floquet states reduces to the straightforward task of diagonalizing 
the rotated Hamiltonian and transforming back to the original frame by applying U + (t). As we shall see later on, the 
natural basis which provides a unified description of the weak- and strong-coupling regimes is nothing but the basis of 
Floquet states associated with a zeroth-order Hamiltonian obtained from the original Hamiltonian (1) by decoupling 
the two doublets, i.e., by taking fii 4 , O23 — ► 0. 
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14 = ElCll + (E4 - LO) (7 44 ((714 + 0"4l) 



(2) 
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II. ANALYTIC MODEL 



We start by performing the time-dependent unitary transformation 
,t r 4 



U(t) 



f 1 

exp i / ^ - (Si + E 2 ) <J U + LJ ((T33 + <7 44 ) - fii2 COS (ojt') (<J 12 + <T 2 1) - ^34 COS (wf') (<T 34 + C7 43 ) 
*/ j ^ 



dfc', (3) 



This transformation, which in particular translates the zero of energy to the point {E\ + E 2 ) /2, enables us to absorb 
the most rapidly oscillating terms of the Hamiltonian (1) and leads, after some lengthy algebra, to the following 
rotated Hamiltonian: 



with 



H' = (A{,/2) {cos [20 (t)] {(722 -<tii) + ian[2(f/(t)] (<r 2 i - <r 12 )} 
+ (A£/2) {cos [2<P"(t)] ((744 - a 33 ) + i sin [20"(t)] (a 43 - <J 34 )} + (A - w) (<7 33 + (744) 
-± (fi 14 + Q 23 ) cos (wi) { e"^ 4 (cos [0_(t)] (<t 23 + <Ti 4 ) - i sin [0_(t)] (ct 13 + cr 24 )) + h.c.} 
•| (f2 23 - rJi 4 )cos (wf) { e _iwt (cos [<j>+(t)] (<7 23 - ctu) - zsin [<f>+{t)] {a 13 - a 24 )) + h.c.} , 



&'(*) = (fi 12 /w) sin (wt) , 



(4) 



(5a) 



b"{t) = (fi 34 /w)sin (ujt). 



(5b) 



0±(i) = ^(t)±0"(t). (5c) 
Next, we express the time-dependent coefficients of H' as a series of Bessel functions J n by using the expansions [36] 

+ OC 



cos [Csin (cot)] = Jo (C) + 2 ^ J 2 „ (C) cos (2nu)t) , 



(6a) 



n=l 



+00 



sin [C sin (cut)} = 2 ^ J 2n+1 (C) sin [(2n + . 



(6b) 



n=0 



This enables us to write the Hamiltonian (4) as a sum of a dominant constant contribution H' and an oscillating 
time-dependent part H[(t), with 



H = ^ Aq R (<t 22 - a n ) + ±A' ,R (cr 44 - cr 33 ) + (A - w) (<7 33 + a 44 ) 

-|^14 ( CT 14 + <?4l) - (^23 + ^32) , 



(7) 



#((*) = A o EiTi ^2n (CO cos (2nwt) (ct 22 - a n ) + iAj, £+^0 -W (CO sin [(2n + l)wt] (<7 2 i - a 12 ) 
+A^' E^i ^2n (CO cos (2nwt) (a 44 - a 33 ) + iA' > -W (CO sin [(2n + l)wt] (<r 43 - ,734) 

I (fi 23 + fii 4 ) [Xc(C-) ( CT 23 + CT14) - %(C-) (°13 + ^24) + h.C.] 



(8) 



I 2 



(Q 23 - Ct 14 ) [Xc(C+) (ct 2 3 - ^14) - %(C+) (0-13 - ^24) + h.c] 



where C = 2Qi 2 /u;, Q" = 2Q, 34 /u and C± = | (C' ± CO- The renormalized splittings Aq R , Aq R and Rabi frequencies 
fi? 4 , 0^ 3 are field-dependent quantities defined as A(, R = A(, J (2fii 2 /w), A[,' R = A^'Jq (2fl 34 /u), Clf 4 = and 



nl\ = ^ R with 



Oi 2 — n 34 \ ui 

and the coefficients Xc(C) an d Xs(C) arc defined by 



qYI _ ^ _J O14 + ^23 j ( ^12 — ^34 \ _j_ ^14 ~ ^23 j /" ^12 + ^34 



Q,\ 2 + il 34 



Xc(C) = ^- 2lut {Jo (C) + e- 2 ^ J 2 (C)} + £ {1 + (1 - M e- 2 ^ 4 } J 2 „ (C) cos (2n^) 



(9) 



(10a) 
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Xs(C) = (1 + e~ 2iut ) J 2«+i (C) sin [(2n + l)wt] , (10b) 

n=0 

where <5„i is the Kronecker delta. 

The important point is that the contribution of the oscillating Hamiltonian H^{t) to the dynamical evolution of 
the system becomes negligible and can be safely neglected under rather general conditions. To see this, we write the 
evolution operator associated with H'(t) = H' Q + H[{t) as the perturbative expansion 

U'(t,0) J dt 1 e iH '^' B' x {t')e- iH ^ +...| . (11) 

It can be easily seen that the integral on the right hand side of (11) is a sum of terms of order A' Q /cu, Aq/cj, 
H14/0J, and 023/^- Thus, for a driving field quasiresonant with the transitions between the lower and upper doublets 
(w « A » Aq, Aq) and weak enough so that the Rabi frequencies of the weak transitions (r2 14 , O23) remain small 
compared to u>, the contribution of H[(t) can be legitimately neglected. This approximation is applicable regardless of 
the value of the coupling constants f2i2 and ^34 and, therefore, is valid in both the perturbative and nonperturbative 
regimes. In particular, in the weak-field regime it leads to the same results as the usual RWA and, consequently, can 
be considered as a nonperturbative generalization of the latter. Under the above conditions, the dynamical evolution 
becomes governed by the Hamiltonian H' which, by defining renormalized energies Ef = — Aq R /2, — Aq R /2, 
Ef ~A — Aq R /2, and i? R = A + Aq R /2, it takes the same form as the weak-field Hamiltonian previously considered. 
Specifically, one obtains H' = H' 14 + H 23 with 

H' m = E?a n + (E? - u) (744 - -^ (<ri4 + <r 4 i) , (12) 

and a similar expression for H' 23 replacing 1 — ► 2 and 4^3. The Schrodinger equation associated with H' can be 
now readily solved analytically to obtain the nonperturbative general solution in the rotating frame 

lv>'W) = £4(i)|i>, (13) 

with probability amplitudes 4(i) given by 

c[(t) = jci (0) cos + A- (4(0)«5 R + c 4 (0)Q R ) sin j e<^'^)\ (14a) 



c' 2 (t) = (c 2 (0) cos ( + i (4(0) 5* + 4(0)O R ) sin ( ^t) ) e <^'^)\ (14b) 




c' 3 (t) = \ c 3 (0) cos (^t) -^(4(0)^ - 4(0)O R )sin(^) \A^ a ~^> i (14c) 



^23 



c' 4 (t) = |4(0) cos ^ftj - ^_ (4(0) 6f 4 4(0)Of 4 ) sin j j e K^~^)\ (14d) 

where we have defined field-dependent renormalized detunings <5 R = i? R — i? R — u) and 5 23 = -E R — E 2 — u, and 

renormalized generalized-Rabi frequencies fi 14 = J (£1 R ) + (<5 R ) and Q 23 = V (^23) + (^23) • The interesting 
point is that while the system dynamics in the strong-field regime is in general rather complicated, when viewed from 
the rotating frame it becomes essentially the same as that of the weak-field regime. The same result holds true in the 
original nonrotating frame by proper choice of the relevant basis. Indeed, by transforming back one obtains 
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and from this general solution one immediately sees that the probability amplitudes associated with the field-dependent 
states \i'(t)) = U + (t)\i) are precisely those given in Eqs. (14a-14d). ft is therefore clear that the renormalized \i'(t)) 
states constitute the natural basis to analyze the time evolution of the system. In fact, when the system dynamics is 
analyzed in terms of such states the nonperturbative effects induced by the strong driving field can be absorbed into 
a redefinition of the relevant energies and Rabi frequencies in such a way that the system evolves obeying the same 
Hamiltonian in both the perturbative and nonperturbative regimes. 
In terms of the original basis, the states take the form 

|l'(i)> = cos<//(t)|l) + i sin <£'(*)] 2), (16a) 
\2'(t)) = isin^'(t)|l) + cos0'(t)|2>, (16b) 
|3'(t)) = e-™*[cos0 // (t)|3>+isin^"(t)|4)] ) (16c) 

|4'(t)> =e- <w *[isin^"(f)|3)+cos^"(i)|4)]. (16d) 

These states constitute a basis of the extended Hilbert space of t-periodic state vectors [37]. In fact, as already 
mentioned, they are the Floquet states associated with the zeroth-order Hamiltonian obtained from the original 
Hamiltonian (1) by decoupling the two doublets, i.e., by taking fii4, SI23 - * 0. This follows from the fact that, in such 
a case, the states {\i)} become the eigenstates of the rotated Hamiltonian H' — H' 14 + H' 2Z (see Eq. (12) and below). 

Note that in the weak-field regime one has <j>'{t), <f>"(t) <C 1 for any t and as a consequence \i'(t)) — > \i), so that the 
renormalized basis becomes indistinguishable from the original one. Similarly, taking into account that Jo (x) —> 1 
and J\ (x) jx — > 1/2 as x ^ it follows that in such regime the renormalized energies and Rabi frequencies approach 
their corresponding bare values, so that, in the weak-field regime, the above formulation simply reduces to the usual 
one. In the strong-field regime, however, the time evolution of the different bare states becomes strongly coupled by 
the driving field and as a consequence it can be rather involved and very different from that occurring in the weak-field 
regime. In contrast, the time evolution of the renormalized states remains always as simple as in the weak-field regime, 
consisting of Rabi oscillations between |l'(t)) <-> |4'(i)} and |2'(i)) <-> |3'(i)). 

III. NUMERICAL RESULTS 

To verify the predictions of the above analytic model, next we perform an exact numerical calculation. We consider 
a quantum particle in a quartic double- well potential driven by an external periodic field of frequency lu (see Fig. 2). 
Since this potential approaches an infinite value at large distances, x — > ±00, it only admits bound eigenstates [38]. 
Consequently, there is no continuum spectrum and such a model is only adequate for describing physical systems at 
energies well below the continuum threshold. 




FIG. 2. Energy diagram of the lowest lying eigenstates of a quartic double-well potential with D = 4. The splittings of the 
doublets have been exaggerated for clarity. 
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Using convenient dimensionless variables, the corresponding Hamiltonian can be cast in the form [39] 




— Ace cos wi. 



(17) 



The dimensionless parameter D determines the barrier height and corresponds, approximately, to the number of 
doublets below the top of the barrier. In the present study we take D = 4. The frequency of the external field has 
been tuned to the transitions between the first and third doublets. Specifically, we have taken ui = (Eq — E\) with 
Ei being the energy levels in the absence of driving field. The dimensionless field intensity A, on the other hand, has 
been chosen to satisfy the strong-coupling condition f2i2/w — 1, where Qi 2 — A(l \x\ 2) (see below). 

To establish a more clear connection with the formalism of the previous Sections, it is convenient to rewrite the 
Hamiltonian (17) in terms of a basis set {\i}} of eigenstates of the quartic oscillator. We have obtained these eigenstates 
numerically, by diagonalization in another truncated basis set {| </?«}} of harmonic oscillator wave functions with a 
conveniently optimized frequency, following the procedure of Ref. [40] . In this way, one gains a complete knowledge 
of the states |z) = X)n( < Pn|*)l9 3 n) by determining the numerical coefficients (ip n \i). 

In terms of the quartic-oscillator eigenstates, the Hamiltonian (17) takes the form 



where cry = and Ei is the energy of state |z) in the absence of the periodic force. Introducing the notation 

flij = X(i \x | j) = X(j \x\ i) for the coupling constants, the connection with the formalism developed previously should 
now be evident. 

Note that, since x is an odd operator and the parity of the quartic-oscillator eigenstate \i) is (— l) i_1 (with i = 
1, 2, 3, . . .), transitions between \i) <-> \j) are allowed only if (i — j) is odd. 

For a sufficiently weak driving field A, any coupling constant f2y will be much smaller than u> and the system will 
evolve in a weak-coupling regime where all of the allowed transitions are weak. Conversely, for a sufficiently intense 
driving field we would have Ojj ^> ui for any i,j and the system would evolve in a nonperturbative strong-coupling 
regime where all of the allowed transitions are strong. However, as mentioned in the Introduction, in between these two 
limiting cases there exists an interesting nonperturbative strong-coupling regime in which the intra-doublet transitions 
become strong while the corresponding inter-doublet transitions remain weak. This is so due to the fact that, in this 
kind of systems, the dipole matrix elements ~ (i \x\ j) between states within a given doublet turn out to be much 
larger than the dipole matrix elements connecting states lying in different doublets. On the other hand, under the 
above conditions, since inter-doublet transitions are weak, contributions coming from off-resonant doublets will be 
negligible, so that one expects the quartic oscillator to behave, to a good approximation, as an effective four-level 
system. Note, finally, that the nonperturbative regime in which we are interested lies within the range of applicability 
of the analytic formalism of Sec. II [see below Eq. (11)]. 

We have solved numerically the time-dependent Schrodinger equation corresponding to the Hamiltonian (18) by 
expanding its solution in the basis set of eigenstates of the quartic oscillator, and have considered as many states in 
the truncated basis sets so as to guarantee well- converged results. Specifically, for the physical parameters considered 
above, the 20 lowest-lying levels of the quartic oscillator have been included, which are more than enough to guarantee 
convergence. As we shall see, under the above conditions, the dynamical evolution of the system can be described, to 
a good approximation, by a four-level model. 

In what follows we shall denote the two states of the upper doublet as |3) and |4), in accordance with the notation 
used in the four-level analytical model developed in the previous Section. Figure 3 shows the time evolution of the 
populations |(i|^(i))| 2 of the bare states |z) (i = 1, 2, 3, 4) for a system prepared in t = in the ground state. The curves 
plotted correspond to the numerical results obtained by solving the Schrodinger equation with the Hamiltonian (17). 
In Fig. 4 we show the corresponding theoretical prediction, obtained from the analytic general solution \ip(t)) given 
by Eq. (15). It is important to note that the numerical results, unlike the analytical ones, include the contribution 
from all of the energy levels (and not only the contribution from the most relevant four levels). In fact, the slight 
discrepancy between Figs. 3 and 4 is due entirely to this circumstance, as demonstrates the fact that both analytical 
and numerical results become indistinguishable when the numerical problem is also restricted to the four most relevant 
levels. 




(18) 
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Figures 3 and 4 show that under the action of the strong external field all of the bare states become highly populated 
and their time evolution couples in such a way that the population dynamics turns out to be quite different from 
the simple Rabi oscillations occurring in the weak-field regime. In contrast, as Fig. 5 reflects, the populations of 
the renormalized states \i'(t)) evolve in time exhibiting the usual Rabi oscillations of the weak-field regime. Solid 
lines in this figure correspond to exact numerical results whereas dashed lines correspond to the analytical results 
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obtained from Eqs. (14a-14d). As before, the small difference between analytical and numerical results originates from 
corrections to the four-level approximation. Indeed, for the high field intensity considered above, the contribution 
of the second and forth doublets to the dynamical evolution of the system, although small, it is not completely 
negligible. In fact, by monitoring the different numerical populations it can be seen that a small proportion of the 
populations of the first and third doublets is rapidly transferred to their corresponding adjacent doublets, giving rise 
to the rapid oscillations that appear superimposed to the usual Rabi oscillations in Fig. 5. When the numerical 
problem is restricted to the four most relevant levels this population transfer vanishes and, as already mentioned, 
both analytical and numerical results become indistinguishable. Since the contribution of level \i) to the dynamical 
evolution of level \j) is proportional to ttij/Sij (where flij is the field-dependent coupling constant between \i) and 
and Sij is the detuning of the corresponding transition [41]), one expects a better agreement between analytical and 
numerical results for smaller field intensities. This is indeed the case as can be appreciated from Fig. 6. This figure 
shows the population dynamics of the rcnormalized states for an external field of the same frequency as before but 
a smaller field intensity, which now satisfies the strong-coupling conditions O12/W = 0.75 (Fig. 6a) and Q12/W = 0.5 
(Fig. 6b). As expected, as the intensity of the external field decreases the four-level approximation becomes more 
and more exact. 




cot 

FIG. 5. Dimcnsionless time evolution of the populations of the renormalized states \i'(t)) for a system initially prepared in 
the ground state. Solid lines are exact numerical results and dashed lines are the analytical results. 

Figures 5 and 6 show that, for the above initial conditions, states |2'(f)) and |3'(£)) remain unpopulated while the 
system population undergoes Rabi oscillations between the renormalized states and |4'(i)), and this occurs in 

both the weak- and strong-field regimes. 
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IV. CONCLUSION 



In the nonperturbative regime the dynamical behavior of driven quantum systems becomes, in general, rather in- 
volved. In this paper we have considered an important class of driven four-level systems which are relevant in the 
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description of numerous processes in molecular and solid-state systems, and we have shown that their nonpcrturbative 
time evolution, when analyzed in terms of a natural basis of renormalized states, essentially reduces to the correspond- 
ing time evolution in the weak-field regime, exhibiting simple Rabi oscillations between the different relevant quantum 
states. 

Such renormalized basis enables one to absorb the nonperturbative effects induced by the strong driving field 
into a redefinition of the relevant energies and Rabi frequencies in such a way that the system evolves obeying the 
same Hamiltonian in the perturbative and nonperturbative regimes. This basis thus provides a unified description 
valid in both the weak- and strong-coupling regimes. In particular, in the weak-field regime, the renormalized basis 
becomes indistinguishable from the original one and the renormalized energies and Rabi frequencies approach their 
corresponding bare values, so that, in this regime, our formulation leads to the same results as the usual RWA and 
thus can be considered as a nonperturbative generalization of the latter. 

This work has been supported by Ministerio de Ciencia y Tecnologi'a and FEDER under Grant No. BFM2001-3343. 
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